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Abstract 

We define a family of string equations with perturbative expansions that admit an 
interpretation as an unoriented minimal string theory with background D-branes and R— R 
fiuxes. The theory also has a weh-defined non-perturbative sector and we expect it to have 
a continuum interpretation as an orientifold projection of the non-critical type OA string 
for c = 0, the (2, 4) model. There is a second perturbative region which is consistent with 
an interpretation in terms of background R-R fluxes. We identify a natural parameter 
in the formulation that we speculate may have an interpretation as characterizing the 
contribution of a new type of background D-brane. There is a non-perturbative map 
to a family of string equations which we expect to be the c = type 0B string. The 
map exchanges D-branes and R-R fluxes. We present the general structure of the string 
equations for the (2, 4&) type 0A models. 
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In what follows, we define a family of theories that have perturbative expansions that admit 
interpretations as unoriented type OA string in the (2, 4/c) superconformal minimal model back- 
grounds, with D-branes and R-R fluxes. The definitions have a well-defined non-perturbative 
regime. For the case of ê = 0, pure supergravity (the (2,4) model) we will see that there is a 
map to a system which has already been found in the literature, arising from studying self-dual 
unitary matrix models in the double scaling limit [1]. We argue that it is natural to identify 
this system as an unorientable type 0B string, the pure supergravity case for that theory. 

The string equation for the k = 2 member of the bosonic minimal string theory for the (2, 2k — 1) 
series can be written as [2-4]: 

-lf" + f = z, (1) 

where a prime denotes vdjdz. This is the Painlevé I equation. Expanding about the large 
positive z regime gives a well defined perturbation theory representing an oriented bosonic 
string theory that is in fact pure gravity, i.e. the Liouville theory, since c = for the matter 
theory. The free energy is related to / as F = —f. This model can be realized as the double 
scaling limit of a one Hermitian matrix model. 

In refs. [5-8] , symmetric matrix models (orthogonal and symplectic) were studied in the double 
scaling limit, and these gave rise to continuum models with contributions from non-orientable 
world sheets. The physics is then encoded in two functions. The oriented contributions come 
from / as before, which (for the k = 2 case) satisfies equation (JTJ), and the unoriented con- 
tributions come from a new function, g, which satisfies (for this k = 2 case) the equation 1 : 

g 3 + Qgg' + 4g"-Qgf-Qf' = 0. (2) 

Given the solution of (JTJ), equation (J2J yields three possible solutions 2 for the asymptotic ex- 
pansion of g(z). We have: 

fM r y2 

J[Z> VZ 24z 2 1152^/ 2 + "' 

/x rx ma v 5u 2 ^/Q 

9lW = ± V6^-_ T _ + ... 

v 2hv 3 15745z/ 
9ÁZ) ~ ~Yz~2A~^~~TE2~~ + '" (3) 

and the free energy of the model is given as the sum of an oriented contribution and an 



^Note that we have taken the liberty of changing the sign of f(z) from the conventions used in [5,7]. This 
is to make contact with the more widespread convention in the literature. 
2 The meanings of which will be discussed later in this paper. 
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unoriented one: 

F = F + F U , where K = ~\í , K = ~9 ■ (4) 

Half of the oriented theory's free energy makes up the oriented contribution to the free energy 
of the model. The remainder comes from g. Unitary matrix models also have critical points, 
and in a double scaling limit performed in ref. [9,10] for the simplest critical point, the following 
equation can be derived: 

h 3 - h" - hz = , (5) 

which is the Painlevé II equation. For large positive z, h has an expansion with an interpretation 
in terms of oriented world sheets. In fact, this model has recently been understood to be pure 
supergravity for the type 0B string theory [11], the k = 1 member of the superconformal series 
(2,4fc). 

In ref. [1] self-dual unitary matrix models were studied. It was found that there are again two 
contributions to the partition function, one coming from orientable surfaces, h, which satisfies 
equation (jSJ; and the other, unorientable contribution, from a function, w, which satisfies the 
Riccati-type equation 3 : 

/ X 9 3 9 / r*\ 

w = -w + z — -h . (6) 
Again, the free energy splits into two parts: 

F = F + F U , where F% = -\h 2 , K=\™ ■ (?) 

It is natural to suppose (and so we propose that it is the case) that this is an unoriented string 
theory based upon some projection of the 0B theory for c = 0. Half of the free energy of the 
oriented theory makes up the oriented sector of this new theory, and the rest is made up of 
unoriented contributions. 

As is by now very familiar, the physics of the even k minimal bosonic string models (of which 
equation (JTJ) is the case k = 2) suffers from a non-perturbative instability, as can be seen from 
the well-known behaviour of solutions of the equation for /. The unoriented theory simply 
inherits this behaviour. 

It is known that there is a family of models that have the same perturbative behaviour as the 
string equations for the bosonic minimal models; but have better non-perturbative behaviour, 
3 We have changed the conventions of ref. [1] by rescaling their t to z/2, their w to w/2, and their h 2 to 2h 2 . 
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defining a family of non-perturbatively stable oriented minimal string theories with the same 
perturbative content. An example equation is: 

fK 2 - ~KK" + l -{Uy = v 2 Y 2 . (8) 

where 

K = -\f" + f 2 -z, (9) 

and for T = they have the same large positive z expansion as equation (JTJ). This follows from 
the fact that K = is a solution to the equation when T = 0. There are solutions that solve the 
equation (JHJ) which are not solutions of K = however, but nevertheless have the same large 
positive z asymptotic expansions as the solutions of K = 0. These solutions define theories 
with well-defined non-perturbative physics [12,13]; and they have now been identified [11] as 
the type 0A minimal string theories, for the (2, 4k) series. The case we have given above is 
in fact k = 2. The string equation above captures more than that however [14]. For r / 
it includes a new feature. For large positive z, the parameter T corresponds to there being T 
D-branes in the background [14], while for large negative z it corresponds to T half-units of 
R-R flux [11]. 

The simplest model is in fact the choice k = 1, which is pure type 0A world-sheet supergravity, 
and has: 

K = f-z, (10) 

In fact, the case k = 1 is very special. There is a change of variables [15] which takes the 
equation (jHl) (with K given in equation (fTTH) . and with T = (for now)) and maps it to the 
equation (JHJ). The map is: 

f = h 2 + z, (11) 

together with the exchange z — > —z. This means, as first pointed out in ref. [11], that the 

physics of pure 0A supergravity and pure 0B supergravity are non-perturbatively related. Non- 

perturbative because these systems each have two distinct asymptotic regimes, large positive 

z and large negative z which are separated by a strongly coupled regime. The map above 

exchanges these regimes 4 . More generally, for non-zero T, the map in equation (JTTj) . together 

with z — > —z, takes the 0A string equation (JHJ) to a generalization [11] of the string equation (JHJ): 

4 As pointed out in ref. [16], this relation may be an important example. The large positive z physics of the 
0A model is identical to that of the k = 1 bosonic minimal string model, which is a topological string theory. 
Pure 0A can therefore be thought of as a non-perturbative completion of the topological theory. Moreover, 
this completion can be interpreted as a 0B model. It is a strong-weak coupling duality for topological strings, 
which deserves further investigation. 
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u 2 Y 2 

h 3 - h" -hz- -^- = , (12) 

which, for the type OB model now, has an interpretation in terms of having F background 
D-branes in one perturbative regime and T units of flux in the opposite perturbative regime. 

A natural question to ask is whether these OA equations have a generalization to include an 
interpretation of the presence of non-orientable surfaces in (at least one of) the perturbative 
regimes. We shall now show that there is a very natural generalization to include this interpre- 
tation. 

First, we digress a short while to introduce an elegant differential operator language in which 
the non-orientable world-sheet physics of the bosonic matrix models can be cast [7] . The basic 
object in terms of which everything follows is the operator: 

d 

Q = d 2 — f(z) , where d = v— . (13) 

For the /cth bosonic model, one forms the following truncation of a fractional power of Q, which 
we shall denote [17] as P: 

P = Q k ~ l • (14) 

The fractional power is defined by requiring that there exist an operator Qï such that it squares 
to give Q. In doing this, one needs a definition of d~ l , and this is such that [18, 19]: 

oo 

d- l f = fd- 1 + ^(-1)7^ , (15) 

i=l 

where means the ith z-derivative of /. The plus sign subscript in the equation (|14j) above 
means that we take only positive powers of d. The string equations are formulated as the 
realization of the fundamental relation [17]: 

[P,Q] = 1, (16) 

where one integrates once with respect to z to get the final string equation. In the case of k = 2 
we have: 

P = d 3 -^{dJ}, (17) 

which leads to equation (JTJ). The unoriented sector is defined in terms of the above objects as 
follows. Require that P be factorizable into the following form: 

P = (d 2 + l -gd+h)(d- l -g) , (18) 
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leading to two equations, which when combined give equation (J2J). We point out here that 
equation (J2J) can be written as 

T' + gT = , (19) 

where 

T = g+\g>- 3 -f. (20) 

A natural sub-family of solutions are those which satisfy T = 0, which is the following Riccati- 
type equation: 

f = \g 2 + \g ■ (21) 

In fact, as pointed out in ref. [5], studying the equation (J2TJ) is equivalent to the more restrictive 
factorization of P: 

P=(d+ l -g)d{d- l -g) . (22) 

It turns out that the gi(z) perturbative expansions from equation (JBJ) are the solutions to the 
T= equation; whereas g 2 (z) corresponds to a solution that satisfies equation (fTÏÏJ). but not 
T= 0. It is actually very interesting that there are these three different solutions for g(z). In 
refs. [5,6], one of the g\ solutions is associated with the orthogonal matrix ensemble; the other 
is associated with the symplectic ensemble. Although the third solution, g 2 , is a solution of 
an equation derived from the matrix model, it has yet to receive an interpretation in terms of 
a particular ensemble. However, there is good cause to expect that it may also correspond to 
viable physics. 

The same operator language that was used above in the [P, Q\ = 1 bosonic case can also be 
used to formulate the type 0A models. This is because, as shown in ref. [12,13], the differential 
operator structure which underlies the minimal bosonic string equations can be used to define 
the minimal type 0A string equations as well. We then apply a factorization condition to the 
appropriate operator involved in this new definition. It works as follows: 

The first derivative of the string equation arises by defining, for a particular k, the operator [12]: 

P = Q + +ï - \zd , (23) 
and imposing the fundamental equation: 

[P,Q]=Q. (24) 

Indeed, this is an equation stating scale invariance, in contrast to the previous one (equa- 
tion (JTEJ)) which states translation invariance. This fits with the fact that the scaled eigenvalue 
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distribution in the underlying matrix model is defined on M + in one case, and R on the other [13]. 
Equation (|24)l then gives the first derivative of the string equation flBJ), and v 2 Y 2 arises as an 
integration constant 5 . 

Now we propose that it is natural to define the unoriented contribution by requiring that P 
factorizes in the same way as P. This will define physics for all the different k cases in a manner 
analogous to the multicritical points of ref. [7], (which we propose to be the type OA (2,4/c) 
unoriented models) but let us work with the case of k = 1 (pure supergravity) to see how things 
work explicitly. 

Before proceeding, we note that the scaling operator defined in equation ()23|) is more restrictive 
than necessary, and misses potentially important physics. More generally, adding a constant, 
vC to the operator P gives an operator which is just as good as a scaling operator, the 
equation ("2*4")) remaining true. 

So we have 6 : 

P = -d 3 + \fd + \f - \zd + vC = -{d 2 + \gd + h){d- \g) , (25) 
and upon expanding and equating powers of derivatives we find two equations: 

-9' + h--/ = -~/ + \z , \hg + \gg' + \g" = \f + vC , (26) 

from which, after elimination of h we obtain: 

g 3 + Qgg' + Ag" - Qgf - Qf' + 2zg = 8uC . (27) 

Given the experience of the previous two cases, the question arises as to whether there is a first 
order Riccati-type equation to which this is equivalent. Noting that ref. [7] showed that one 
can restrict to that form by a more specific factorization of P, we try the same thing for P: 

P = -d 3 + ~fd + lf - \zd + vC = -{d + \g)d{d - \g) , (28) 

giving: 

- -/ -g' = -\f + \z , \gg' + \g" = \f + uC (29) 

In order for these two equations to make sense we require that they must have the same content. 
To do this we must set C = —1/4, which means the equation for g can be written simply as: 

S = , (30) 

5 There are many other ways to derive equation {HJl. See refs. [12, 13] for matrix model derivations 
6 Note that we have chosen the negative sign of the square root of Q here to make things appear more natural 
later. 



7 



where 

S = g +\g 2 -\f + \z = T+\z. (31) 

Using this notation we find that the result of the less restrictive factorization ()27|) can be 
expressed as (again with C = —1/4): 

S' + g S = , (32) 

This is very pleasant indeed, since using the map together with z — > —z and the iden- 
tification g = w, the equation S = becomes equation (jHJ) derived for the symmetric unitary 
matrix models, which as we said earlier pertain to the type 0B system! So we have extended 
the map between 0A and 0B to one between their unoriented descendants. This also serves to 
put our computations on a firm double scaled matrix model footing. 

It is interesting to examine the behaviour of the solutions in each asymptotic regime. We have 
solutions of equation (jHl) for f(z) in the large positive z and negative z directions, respectively: 

5t/ 3 r(4r 2 + i) z/ 4 r 2 (8r 2 + 7) 

U[Z) - Z + ~z~^ + 32z^ 8^ + "' 



= , 2 ( 4 r 2 i) / + 4 , 2 (4r2 - 9) + ^ (4 r 2 - 9) (4 r 2 - 21) + \ 

Az 2 \ z 6 lQz b ) 

Looking in the positive z expansion, for example, integrating twice and dividing by v 2 reveals an 
expansion in the dimensionless parameter g s = u/z 3 ^ 2 . The first term apparently corresponds 
to a sphere term, but it is non-universal and should be dropped, while the next two correspond 
to the contributions of the disc and annulus, respectively. The relevant diagrams come with a 
factor g 2h+b ~ 2 T b , where b is the number of boundaries and h is the number of handles. In the 
negative z expansion, T 2 is to be taken as representing a single R-R flux insertion [11], and so 
the first term is the torus, (the T independent term) and sphere with a flux insertion. 

There are unique non-perturbative completions of these solutions, and some of these are plotted 
in figureHfor a range of values of T. With each of these asymptotics comes two choices for g(z), 
as solutions of S = 0. For the large positive and negative z directions, we have the choices, 
respectively: 

, nr . u(±3T-2) í/ 2 (±2ir 2 -36r±20) 
9+ , l( z) = ±2^ + L_ 1 .^ _ > + ... 

, . r- r- v iV2v 2 (12T 2 - 13) 3z/ 3 (12T 2 - 13) 
9-M = ±«^--=F ^ 1 4^ L + --- (34) 
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f(z) 




Figure 1: The non-perturbative solution for the function, f(z), which contributes to the oriented 
sector of the model, for a range of positive values of T. The bottom curve is the T = case. 



and we can also have third and fourth solutions that satisfy (jH2j) but not (JÏÏUJ) : 



9+,2(Z) 



v 9z/ 2 r 15z/ 3 (3r 2 + 2) 

+ ~ ~7~ ~~~1 " + ■ 



9-,2( Z ) 



z 

V 

z 



4 z 5/2 8z 4 

3z/ 3 (12 T 2 - 13) 3z/ 5 (336 T 4 - 3384 T 2 + 3085) 



16*' + - (35) 

An examination of the positive z regime for the first case (g+.i) shows (after integrating once 
and dividing by v) an expansion in g s again. It is clear that the first two terms are the 
contributions of RP 2 , the Mobius strip, and the Klein bottle, respectively. The diagrams come 
with a factor g 2h + b + c - 2 r b where c is the number of crosscaps. In the negative z regime we have 
an interpretation with just R-R fluxes again, the first term being MP 2 and the next the Klein 
bottle. The next has contributions from closed surfaces with the following: c = 1, h = and a 
flux insertion; c = 3, h = 0; c = 1, h = 1. 

For the second case, there is again a surface interpretation consistent with there being D-branes 
in the positive z regime, and fluxes in the negative z regime. However, there are missing orders 
that deserve some explanation, which is perhaps to be found in a study of the continuum 
theory 7 . 

7 There have been some recent studies of non-critical string theory on unoriented surfaces which may be 
relevant [20-23]. 
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Plots of the non-perturbative completions of these two cases are given in figure 121 As in the 
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Figure 2: Two non-perturbative solutions for the function g(z) which contributes to the unoriented 
sector of the model, for a range of positive values of V. In both cases the top curve is the T = case. 

case of the [P, Q] — 1 formulation, it is tempting to speculate that all of these perturbative 
expansions for g(z) are on equal footing; that is, they are equally valid in non-orientable string 
theory Whether or not all these models stem from matrix models is also a question worth 
investigating. 

Having arrived at the case Co = —1/4, it is natural to wonder whether the other choices for 
Co are physical as well. After a little algebra, we find that we can rewrite equation (}2*Tj) as the 
following: 



where we have defined the constant G for later convenience. The case before was G — 0. To 
what do other values of G correspond? The case of Cq = (G = 1/2), which corresponds to 
the P operator originally defined without the shift, gives the following expansions for g (again 
two choices in each regime). For large positive z we have: 




(36) 



9+,Á z ) 




2 (±84T 2 - 120T±47) 

QAz 5 / 2 



+ ••• 



9+Á z ) 



3v 3z/ 2 r 3z/ 3 (76 T 2 + 77) 



_l_ . . . 



(37) 



2z z 5 / 2 32z A 
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while for the large negative z regime we have: 

r r- 3v %Á2v 2 (24 T 2 - 47) 21z/ 3 (8 Y 2 — 17) 

g _ M = J^-l)_6^(8r'-MI« + 13) + .„ (3g) 

Again, a surface interpretation can be given along the lines of those discussed below equa- 
tions (jHl) and (jSSI). 

Leaving G as a parameter in the solutions (rather analogous to T, which controls the number 
of background D-branes or R-R fluxes), these equations yield the following solutions for g[x)\ 

r v(±3T-2 + G) z/ 2 (±2ir 2 + 12r(G-3)±3G 2 + 18G±20) 
g +tl (z) - ±2Vz+ ^ + •■■ 

r- u(G + l) iVïv 2 ( 12 T 2 -QG 2 - 18 G- 13) 
9-Á*) = W2VS- K — ' + ^ '- + ■■■ 

v(G + l) 3z/ 2 r(2G + 3) z/ 3 (3r 2 (8G+ 15) + 2G 2 (G + 9) +46G + 30) 

9+ÁZ) ~ z + 4i^ 8? + "' 

z/(2G-l) z/ 3 (24r 2 G-36r 2 -16G 3 + 72G' 2 -98G' + 39) 

9-á z ) - - + ^4 + " ' y ) 

These solutions interpolate between those displayed earlier for the cases G = and G = 1/2. 
Figure El shows some sample plots of the full non~perturbative g(z) for varying values of G. 
Upon examination of the expansions, it is tempting to interpret G as a parameter controlling 
some aspect of the background, which is present only for the unoriented case. One possibility is 
that it incorporates the presence of some number, G, of a new type of D-brane that is present 
only in non-orientable diagrams (since it does not appear in the contributions from oriented 
surfaces). A candidate such brane could be one which is forced to remain stuck at (the analogue 
of) orientifold planes, a sort of fractional brane. This deserves further investigation. It is the 
case G = that has a known matrix model interpretation (by virtue of the map to self-dual 
unitary matrix models discussed earlier), and this would correspond to having such branes 
absent. 

A possible clue to the nature of G might be found in an examination of the instanton corrections 
to the expansions (terms exponentially small for small g s ), which often have an interpretation 
as D-branes in the theory [24,25]. Some standard computations reveal that there are a number 
of instantons associated with each solution. These have actions (D-brane tensions) given by: 

41 41 21 21 21 

f+ ■ o — 5 9+,i ■ o— or ~q ' 9+,2 ■ -— ; g+ ■ -— ■ (40) 
3g s 3g s 3 g s 3 g s 3 g s 
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Figure 3: The non-perturbative solution for the function g(z) which contributes to the unoriented 
sector of the model, for a range of values of G. The middle curve is the case of G = 0. 

In each case, the positive z instanton has a negative z counterpart which has an action precisely 
1 / y/2 times smaller 8 . Here, the solution g + is the one obtained directly from the equation S — 0. 
In each perturbative regime it appears that the tension of the unoriented expansion instantons 
takes one of two values: the same value associated to the oriented instanton; and half the 
value associated to the oriented instanton. It is therefore possible that this second unoriented 
instanton is the extra D-brane that is in the background, whose number is measured by G. This 
needs further work to properly establish however; but our picture is suggestive, if somewhat 
speculative. 

The functions which give rise to / and g in the matrix model arise from orthogonal polynomial 
coefficients that come from effectively splitting the eigenvalues Aj of the matrices into the 
ones with even labels (\2j) and the ones with odd labels (\2j+i), respectively. The physics 
that arises from quantities controlled by / are known to have an interpretation in terms of 
an eigenvalue distribution, the properties of which in the scaling limit can be interpreted in 
terms of background D-branes. In fact there is a direct map between the function / and the 
eigenvalue distribution which can be studied on the sphere using an integral representation [26], 
or in terms of the resolvent of the operator Q [27]. The latter is particularly useful for going 
8 It is not clear if these latter represent branes given the R-R flux interpretation in that perturbative regime. 
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beyond the sphere and incorporating the effects of T in the background, as recently shown in 
ref. [16]. It is natural to wonder if there is an effective eigenvalue distribution associated to 
the function g that controls the physics of the unoriented sector. The instantons we computed 
above would have a natural home there in terms of the value of an effective potential at an 
extremum, and may correspond to a new type of D-brane. It seems likely (since g arises as 
part of the operator P which is made of a fractional power of Q) that there is some natural 
geometrical information to be extracted from g analogous to an eigenvalue distribution and its 
effective potential. This is worth exploring. 

Finally, it is known [16,28] that the introduction of T D-branes (or R R flux units) into a closed 
minimal string background can be performed very elegantly in the language of the underlying 
twisted boson (</?, which is related to the loop operator and the resolvent) by simply acting 
with a vertex operator V =: exp(— r <p/V2) : It would be interesting to see how the unoriented 
theory is described in this framework. Perhaps the role of G may be best uncovered in this 
context. 

Acknowlegements: JEC is supported by an EPSRC studentship, and also thanks the Uni- 
versity of Southern California's Department of Physics and Astronomy for hospitality. 
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